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We use the moduli space approximation to study the time evolution of mag- 
■ netically charged configurations in a theory with an SU(iV + 2) gauge symme- 

try spontaneously broken to U(l) xSU(iV) xU(l). We focus on configurations 
, containing two massive and N — 1 massless monopoles. The latter do not 

appear as distinct objects, but instead coalesce into a cloud of non-Abelian 
^| field. We find that at large times the cloud and the massless particles are 

i-^H ■ decoupled, with separately conserved energies. The interaction between them 

J> ■ occurs through a scattering process in which the cloud, acting very much like 

a thin shell, contracts and eventually bounces off the cores of the massive 
monopoles. The strength of the interaction, as measured, e.g., by the amount 
of energy transfer, tends to be greatest if the shell is small at the time that 
it overlaps the massive cores. We also discuss the corresponding behavior for 
the case of the SU(3) multimonopole solutions studied by Dancer. 
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I. INTRODUCTION 



There has long been an interest in the magnetic monopoles that arise as classical solutions 
in certain spontaneously broken gauge theories. When the theory is quantized, these give rise 
to magnetically charged particles that can be regarded as the counterparts of the electrically 
charged elementary quanta of the theory. In certain supersymmetric theories, these two 
classes of particles are believed to be related by an exact duality symmetry [[]]]. 

When the gauge symmetry is spontaneously broken, the gauge bosons corresponding to 
the generators of the unbroken subgroup remain massless. If the unbroken subgroup contains 
a non-Abelian factor, some of these massless gauge bosons transform nontrivially under the 
gauge group and thus carry an "electric" charge. Duality suggests that these should have 
massless magnetically-charged counterparts. Although these cannot be realized as isolated 
classical solutions, evidence for their existence can be found by analyzing multimonopole 
solutions. In this paper we will investigate some of the properties of these massless monopoles 
by examining the role that they play in low energy scattering processes.^ 

Recall that an arbitrary magnetically charged Bogomol'nyi-Prasad-Sommerfield (BPS) 
@ solution can be naturally understood as being composed of a number of fundamental 
monopoles of various types ||. If the gauge symmetry G is broken to an Abelian subgroup 
U(l) r , as happens for generic choices of the adjoint Higgs vacuum expectation value, there 
is an integer-valued topological charge for each U(l) factor. Associated to each of these is 
a fundamental monopole carrying a single unit of that topological charge. Each of these 
fundamental monopoles can be realized as a classical solution by embedding the unit SU(2) 
monopole solution in an appropriate subgroup of G. 

For the special choices of the Higgs vacuum expectation value that give a non-Abelian 
unbroken subgroup, the BPS mass formula implies that some of the fundamental monopoles 
should become massless. Although the corresponding embedding solutions reduce to pure 
vacuum solutions in this limit, analysis of multimonopole solutions suggests that the col- 
lective coordinates of the massless fundamental monopoles survive as degrees of freedom. 
Examining static classical solutions with magnetic charges corresponding to a sum of mas- 
sive and massless monopoles, one finds that the massive monopoles are surrounded by one 
or more "clouds" of non-Abelian field ||. The number of collective coordinates || needed 
to describe these clouds is exactly equal to the number that one would have attributed to 



Although the underlying motivation for this work arises from theories with extended supersym- 
metry, this supersymmetry does not come into play at the leading, classical, order to which we 
work. 
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the massless monopoles.0 



To make these ideas more concrete, consider the case of an SU(iV + 2) gauge theory with 
an adjoint Higgs field whose asymptotic value (in some fixed direction) can be brought into 
the form 



with t\>t2>... tjv + 2- The asymptotic magnetic field in the same direction can be written 
as Fy = eijkTkQM/r 3 , where the magnetic charge Qm is of the form 



The rik are all integers, and give the number of fundamental monopoles of the various types. 
The mass of the kth species of fundamental is proportional to tk — tfc+i- Hence, if the tj are 
all distinct, so that the symmetry is broken maximally, to U(1) 7V+1 , the N + 1 species of 
fundamental monopoles are all massive. If M > 2 of the tk are equal, there is an enlarged 
unbroken symmetry, with M — 1 of the U(l) factors being replaced by an SU(M). 

We will focus on solutions containing one of each species of monopole, with t\ ^ — 
£3 = . . . = t jv+i 7^ tN+2- The unbroken group is then U(l) xSU(iV) xU(l), and our solutions 
contain two massive and N — 1 massless monopoles. This example has two advantages for 
us. First, a number of analytic results about the static solutions and their moduli space 
are already known. Second, because the two massive monopoles correspond to commuting 
SU(2) subgroups of the SU(A" + 2), they can have no direct interactions. Any interaction 
between them must be due to the mediation of the massless monopoles, and thus can give 
us insight into the properties of the cloud. 

Explicit solutions for an SU(iV + 2) theory with the Uk all equal to unity were obtained 
in Ref. @. For either type of symmetry breaking, these solutions are described by 4(iV + 1) 
collective coordinates. When the fundamental monopoles are all massive, these have a 
natural interpretation as the spatial coordinates, x a , and U(l) phases of the N+l constituent 
monopoles; as long as the separation between the monopole positions is large compared to the 
monopole core sizes, the multimonopole nature of the configuration is evident in the classical 
solution. As the Higgs field approaches the value corresponding to the non-Abelian breaking, 



2 This counting of collective coordinates assumes that the total magnetic charge of the configuration 
is Abelian. Configurations whose total magnetic charges have non-Abelian components can be 
regarded as having clouds that have expanded out to spatial infinity, and have a corresponding 
reduction in the number of collective coordinates. For discussions of some of the pathologies that 
arise in the presence of non-Abelian magnetic charges, see ||. 



$ = diag(ti,* 2> ---,*!V+2) 




Qm = diag (n x , n 2 - ni, . . . , n N+1 - n N , -njv+i) • 



(1.2) 
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monopoles 1 and N + 1, which remain massive, clearly retain their identity.^ However, the 
N — 1 monopoles that become massless are replaced by a single cloud surrounding the 
massive monopoles; this cloud can be viewed as being formed by the merger of the cores of 
the massless monopoles. Inside this cloud, the magnetic field is approximately equal to the 
the Coulomb field appropriate to the charges of the massive monopoles; it has both Abelian 
and non- Abelian components. Outside the cloud, the Coulomb component of the magnetic 
field is purely Abelian, with the non-Abelian contributions falling with a higher power of 
distance. 

The size of this cloud is characterized by a "cloud parameter" b that is equal to the 
sum of the separations |x a — x a+ i| between the monopoles. Any value b > r = |xi — 
Xjv+i|, the separation of the two massive monopoles, is possible, and the static energy is 
independent of b. All other trace of the massless monopole positions appears to be lost in 
this limit. In particular, the fields show no special behavior at the x a corresponding to the 
massless monopoles. In fact, any change in these "positions" that leaves b invariant can be 
compensated by a gauge transformation of the solution. 

Our goal in this paper is to gain further insight into the nature of these massless 
monopoles by examining the role that the cloud plays in monopole scattering. To do this 
we use the moduli space approximation M , which reduces the dynamics of the infinite num- 
ber of field degrees of freedom to that of a finite number of collective coordinates g^. The 
Lagrangian for this reduced set of variables can be written as 



is a naturally defined metric on the moduli space of solutions. The evolution of the collective 
coordinates is simply geodesic motion in this metric. When the monopoles are all massive, 
this reduction to a finite number of degrees of freedom fits quite well with the point of 
view that particles arising from classical soliton solutions should be seen as having a similar 
status to those arising as quanta of the elementary fields, and that the solutions with higher 
charge should be interpreted as multiparticle states. When some of the monopoles become 
massless, the classical solutions lose their manifest multiparticle nature. Nevertheless, the 



3 We label the component monopoles by the particular magnetic charge that they carry; thus 
monopole a has n a = 1 and corresponds to an SU(2) embedding in the a and a + 1 rows and 
columns of the SXJ(N + 2) matrices. 




where 
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massless monopole collective coordinates survive (although with their physical interpreta- 
tions modified) and the moduli space Lagrangian has a smooth limit J9|, provided that the 
total magnetic charge is Abelian. 

The moduli space approximation is usually expected to be valid in the limit of small 
monopole velocities. This ensures that the monopole kinetic energy is small enough that the 
excitation of the nonzero-frequency modes of the massive fields is energetically suppressed. 
Although simple energetic arguments are not sufficient to rule out excitation of any massless 
fields that might be present, more detailed analysis flO| shows that the approximation holds if 
the long-range fields are all Abelian. As we will see, the situation with massless non-Abelian 
fields is more complex. 

To use the moduli space approximation, we need the metric on the space of multi- 
monopole solutions. The exact form of this metric is known for only a relatively few cases. 
However, these include [0,0] the case where there is at most one of each species of funda- 
mental monopole, which is precisely what we need. Unfortunately, the form of this metric 
that was obtained in Ref. |12| was expressed in terms of the positions and U(l) phases of the 
(N + 1) individual monopoles. Although these are the natural coordinates to use when the 
monopoles are all massive, they are much less useful in the limit we are interested in, where 
the massless monopole positions have little direct physical meaning. A more natural set of 
coordinates would include the positions and U(l) phases of the two massive monopoles, the 
cloud parameter b, and the parameters needed to specify the global SU(iV) orientation of 
the cloud. 

Thus, our first task is to change to a new set of variables. We do this in two steps. 
In Sec. II we transform to an intermediate set of coordinates that allows us to separate 
the metric into two parts, one that gives the terms in the Lagrangian containing r and 
b, and an "angular" part that contains the entire dependence on the spatial and gauge 
orientation angles and phases. The natural next step would be to relate these intermediate 
coordinates to the spatial Euler angles and their gauge counterparts. This turns out to 
be somewhat nontrivial. However, we can bypass this step by rephrasing the problem in 
terms of the corresponding angular momenta and charges. We do this in Sec. Ill, leaving 
the discussion of the Euler angles to the Appendix. With these results in hand, we are 
then ready to determine the geodesies of the metric and the behavior of the monopoles in 
scattering processes. This is described in Sec. IV. We summarize our results and make some 
concluding remarks in Sec. V. 
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II. THE MODULI SPACE METRIC 



The moduli space of solutions for a system comprising two distinct fundamental 



monopoles was obtained in Ref. [II]; the extension to any number of distinct monopoles 
was conjectured in Ref. [T2[] and proven in Ref. Locally, this space can be written as a 
product of a flat 4-dimensional space, corresponding to the center-of-mass position and an 
overall U(l) variable, and a relative moduli space. For the case of interest to us, an SU(iV+2) 
theory with one monopole of each species, the moduli space is 4(iV + l)-dimensional, and 
the 4iV coordinates for the relative moduli space can be chosen to be the relative positions 

= — x^+i (A = 1, 2, . . . , N) and the corresponding relative U(l) phases if) a. 

In Ref. H it was shown that when the SU(N + 2) symmetry is broken to 
U(l) xSU(iV) xU(l) this metric can be written as 



/re! 



-dr 2 A + r A^dtpA + cos 6*^ 



- A 



J2 r A (#a + cos A d<j> A 

. A 

(2.1) 



where va, a , and <pA are the spherical coordinates of the vector r^, /x is the reduced mass 
of the two massive monopoles, and 

,2 
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in 



A 



9 2 V 



2b + r c 



(2.2) 



g 2 + 87r/x& 

Here g = An/e (with e the gauge coupling of the theory) is the magnitude of the fundamental 
monopole magnetic charge, while r c = 2k/ fi is approximately equal to the sum of the core 
radii of the two massive monopoles. 

The difficulty with using this expression for our purposes is that, as we noted in Sec. I, 
the positions of the massless monopoles do not have any direct physical meaning. We would 
therefore like to transform to a more physically meaningful set of variables. This set should 
include the cloud parameter b and the spherical coordinates r, 9, and of the vector r 
that gives the relative separation of the two massive monopoles, but must also include the 
relative U(l) phase and a number of SU(iV) orientation angles. 

One might expect that iV 2 — 1 = dim [SU(iV)] such angles would be needed to specify the 
gauge orientation of the solution. For N > 2 this would require more collective coordinates 
than are available. However, examination of the explicit solutions shows that there is always 
a U(iV — 2) subgroup of the unbroken SU(iV) that leaves the fields invariant [i.e., any of 
these SU(iV + 2) solutions is essentially an embedding of an SU(4) solution]. Hence, there 
are only 4N — 5 = dim [SU(A^)/U(A^ — 2)] global SU(iV) parameters. Together with the 
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U(l) phases, the cloud size b, and the positions of the two massive monopoles, this gives a 
total of 4(iV + 1) collective coordinates, as required. 
We begin by defining complex coordinates^ 

z \ = v ^cos(^ A /2)e- i{ ^ + ^ )/2 

z \ = > /rJsin(0 j4 /2)e i ^-^)/ 2 . (2.3) 

These satisfy 

r A = z\{a) ah z\ (2.4) 

where z a A = (z A )* and the are the Pauli matrices. Hence, the relative position of the two 
massive monopoles is^ 

r = zcrz (2.5) 

while 

b = zz. (2.6) 

When rewritten in terms of these coordinates, the metric takes the form 

Grei — fi [z a- dz + dz a- z] 2 + 4k dz dz + A [z dz — dz z} 2 

= fxdr 2 + Andzdz + X[zdz — dz z] 2 . (2.7) 

Next, we find an SU(2) matrix U satisfying 

Uv ■ aU- 1 = a 3 (2.8) 

where r = r/r. (Thus, U corresponds to a rigid rotation of the entire assembly of monopoles 
that puts the two massive monopoles on the z-axis.) This allows us to define a new set of 
complex variables 

w a A = U ab z b A (2.9) 
in terms of which r and b are given by the simple expressions 

wajW = r5js , ww = b . (2-10) 



4 In the notation of Ref. I], z\ = £a/2 and z\ = Q/2. 

5 Omitted indices on the z\ and related quantities should be understood as being summed over. 
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These identities imply that the w A must be of the form 



1 b + r (!) 

w a = y—2-PA 



where the p A are a pair of complex iV-component vectors obeying 

>U b ) 



PW = 8* ■ (2.12) 

For later use, we also introduce N — 2 orthonormal vectors (r = 3, 4, . . . , N) orthogonal 
to the p A * ■ 

We must now rewrite the metric in terms of these new variables. It will be helpful to 
introduce some new notation. First, we define two unit vectors ni and n2 orthogonal to r 
by requiring that 

Uh a ■ aU- 1 = a a . (2.13) 
Next, we define Aj and Vj (j = 1, 2, 3) by 

dUU- 1 = -—A j a j (2.14) 

and 

Vj — i[w Oj dw — dw Oj w) . (2-15) 

Finally, we define 

Vq = i [w dw — dw w] . (2-16) 
We also need two identities. First, by differentiating both sides of Eq. ( [2.8|) , one is led to 

dr = rdr + r(h 1 A 2 - n 2 ^i) • (2.17) 
Next, viewing the w a as vectors in an iV-dimensional complex vector space, we write 

dw dw = dw a A [Yi L AB + Yi T AB ] dw a B (2. 18) 

where 

nS B = Ep& , pS ) (2.19) 

0=1 



is the projection operator onto the two-dimensional subspace spanned by w 1 and w 2 (or by 
z l and z 2 ) and 



n 



AB 



N 

E 

r=3 



(2.20) 



is the projection operator onto the orthogonal (N — 2)-dimensional subspace. Using 
Eqs. ( |2.10| ), ( p.llj ), ( P-15| ), and (|2.16|) , we obtain the second identity, 



, , \d(b + r)] 2 \d(b-r)} 2 b 
dw dw = ^hz t- + ^ ~~~ — + 



i(b + r) 



EK 2 - 



(b-r) A(b 2 - r 2 ) ^ " 2(b 2 - r 2 ) 



:2 AV V 3 ) + dwU T dw . 



(2.21) 



We have not yet obtained explicit expressions for 9, 4>, or the gauge orientation phases 
in terms of the w\ or U. However, it is easy to see that neither the and V v nor the 
combination dw I1 T dw can contain any factors of dr or db. With this in mind, and using our 
two identities and the above definitions, we substitute Eq. (|2.9| ) into Eq. (|2.7|) and obtain 

12 1 



£rel = fJ>dr + K 



2 j [d(b + r)] 2 [d{b-r)]< 



\ 2(b + r) 2(6 -r) 



+ G: 



ang 



(2.22) 



where 

2 

Gang = E 



i=l L 



Kb 

^2 ty 2 



Kb 



^2 ty 2 



V 2 + k {2A i V l + 



k6 



^2 ty 2 



A (Uo + rAs)^ 



(V 3 + bA 3 ) 2 



2nr 

^2 ty 2 



(K) + rA 3 )(V r 3 + 6A 3 ) +4^«jn T dw (2.23) 



contains the entire dependence of the metric on the spatial angles and the various gauge 
orientation phases. In the next section we will rewrite Q ang in a form that will allow us to 
analyze the evolution of these angle and gauge variables. 

Before doing so, we comment on a point that we have skipped over. Equation ( |2.8[ ) does 
not completely determine the matrix U, but instead leaves the freedom to make a further 
transformation of the form 



U — ► e ixa ' i/2 U . 



This implies a corresponding transformation 



w 



e iX*3/2 w 



(2.24) 



(2.25) 



of the w% and rotates the unit vectors fii and fi2 so that 



n! + m 2 — >e x {n 1 + zn 2 



(2.26) 
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This ambiguity in the definition of U is not simply due to our failure to impose a suf- 
ficient number of conditions. Instead, it reflects the fact that these solutions possess an 
axial symmetry. They are invariant under rotations about the axis joining the two massive 
monopoles, provided that these are combined with an appropriate global gauge transforma- 
tion. Because of this symmetry, there is no natural way to pick out any particular pair of 
axes in the plane perpendicular to f . As a result, the form of the metric should be invariant 
under such redefinitions of U. To verify that it is, we note that the auxiliary quantities Aj 
and V u transform as 

Ax + %A 2 — > e~ ix {A x + %A 2 ) 

A 3 — > A 3 + d X 
V l + iV 2 ^e- ix (V 1 + iV 2 ) 
Vs^Vs-bdx 

Vo-^Vo-rd X . (2.27) 



Substituting these results into Eq. ( |2.23| ), we see that the metric is indeed unchanged, thus 



providing a useful consistency check on our calculations. 

III. ANGLES, PHASES, AND CONSERVED QUANTITIES 



In this section we will bring Eq. Q2.23 ) for (? ang into a form more suitable for analyzing 



the dynamics of the angle and phase variables. The most straightforward method for doing 
this would begin by defining generalized Euler angles to describe the spatial and the U(iV) 
phase orientations. Given these, one can construct a set of 1-forms cu- 7 invariant under both 
the spatial SO (3) and the internal SU(iV) symmetries of the theory, and write 

£an g = Iij(r,b)u l u J . (3.1) 

The symmetries of the system imply the existence of a conserved angular momentum and 
conserved U(iV) charges. The "body- frame" components of these quantities, defined by 

X i dt = I i j{r,b)J (3.2) 

are the variables conjugate to the invariant one-forms. Writing 

G aBS = (r 1 ) ii (r,b)X i X j dt i (3.3) 



and then substituting this into Eq. ( |2.22| ) gives a convenient starting point for studying the 
dynamics of the system. 

A complication in our case is that, because of the extra symmetries of these multi- 
monopole solutions, the number of angular variables is less than might be expected. For 
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most choices of the magnetic charges, an arbitrary configuration would require 3 spatial 
Euler angles and iV 2 phases specifying the U(iV) orientation. By contrast, the relative mod- 
uli space in our case is 4iV-dimensional so that, after extracting r and b, we are left with 
only 4N — 2 angle and phase variables. The "missing" (N — 2) 2 + 1 variables are explained 
by extra invariances of the solutions we are considering. First, we noted previously, these 
SU(iV + 2) solutions are essentially equivalent to embeddings of SU(4) solutions, so that 
there is always a U(iV — 2) subgroup that leaves the configuration invariant.^ The (N — 2) 2 
"Euler angles" corresponding to this subgroup play no role in the dynamics, and hence do 
not enter the Lagrangian. Second, the axial symmetry discussed at the end of the previous 
section implies that some linear combination of a spatial Euler angle and a U(iV) phase must 
also be absent from the metric. [] 

The fact that these extra symmetries involve mixtures of spatial rotations and gauge 
transformations complicates the explicit construction of the uo l appearing in Eq. ( |3.1| ); they 
are not simply subsets of the standard invariant one-forms for SO (3) and U(iV). Fortunately, 
it turns out that we can bypass the construction of the uj % and go directly from Eq. (|2.23|) 
to an expression of the form given in Eq. ( |3.3| ). We will follow this procedure in this section, 
leaving the discussion of the generalized Euler angles and the of to the Appendix. 

We begin by identifying the conserved quantities of the system. Recall that if a metric 
has a Killing vector of the form 

then the charge 

Q (i) = Kfaq a = g ab K a (i) q b (3.5) 

is a constant along any geodesic. Since the classical solutions of the moduli space Lagrangian 
are all geodesies, these Qu\ are the desired conserved charges. 



6 Note, however, that the particular U(iV — 2) subgroup depends on the SU(iV) orientation of the 
configuration, and thus changes as the various phases vary with time. 

7 The apparent analogy between a rigid body symmetric top and our axially symmetric field 
configurations is a bit misleading. In the former case, time-dependent rotations about the symmetry 
axis are physically meaningful and have kinetic energy associated with them. Although the rotation 
angle about this axis does not appear in the Lagrangian, its time derivative does. In our case, 
such rotations, whether time-dependent or not, are undetectable and have no associated kinetic 
energy. Thus, the best analogue is not an ordinary symmetric top with principle moments of inertia 
1% = I2 7^ I3 7^ 0, but rather an infinitely thin top with I\ = I2 and 1$ = 0. 
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The Killing vectors of the metric of Eq. ( |2.1|) were determined in Ref. ||. Written in 
terms of our complex variables z\, the rotational Killing vectors take the form 



Lk 



while the N 2 Killing vectors 



E 



b ( \ 3 -b ( \ ® 
Z A \ a k)ab - Z A [a k ) ba — 



(3.6) 



AB 



V2 



■7°- - — 7 a 

Za dz% Zb 



d 



(3.7) 



satisfy the algebra of U(iV). 

Using Eqs.( p.4|) and ( |3.5| ), we construct from these the angular momentum 



J dt = fi (z a z) x[z cr dz + dz cr z\ + %k\z a dz — dz a z\ — iX [z dz — dz z] (z cr z) 

= fj, r xdr + k i(w UaU' 1 dw - dw UaU' 1 w) + iw(A U cr U' 1 + UaU' 1 A)w 
-Ar (Vo + rA 3 ) 

(where A = Ajtjj) and the U(iV) charges 

Sab dt = iy/2n [dzs za — zb dzA] + iX zb za [z dz — dz z] 

= V2k [i(du>B Wa — wb dwA — wb ctj wa Aj] + a/2A % wa(V + rA 3 ) . 



(3.8) 



(3.9) 



The components implied by the bold-face notation in Eq. (|3.8|) and by the subscripts 
in Eq. ( |3.9|) are those corresponding to a fixed "space-frame" and are conserved. What we 
actually need are the components corresponding to a moving "body-frame" . For the angular 
momentum, this requires identifying three axes that rotate with the monopole system; a 
suitable choice is hi, n 2 , and h s = r. The body frame components (which are themselves 
rotational scalars) are then 



J; dt = n, • J dt 



(fir 2 + nb)Ai + KVi 

-rX(V + rA 3 ) + k(V 3 + bA 3 ) 



1,2 



(3.10) 



In obtaining this, we have used Eqs. ( |2.13| ) and fl2.15| ). 

To construct a "body- frame" for the U(iV) charges, we need a set of iV complex basis 
vectors in the internal space that transform under the fundamental representation under the 
action of the Sab- We can take these to be the and the that were introduced in 
Sec. 2. Using the p£ and the Pauli matrices, we can construct three U(iV)-invariant charges 



V2 



PA [ T k\abPB t-AB ■ 



(6) 



(3.11) 
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(The definition of these has been chosen so that they correspond to generators of an SU(2) 
subgroup with the standard normalization.) Straightforward calculations give 



/v 



Ttdt 



f Vi + Vb 2 - r 2 Ai 



i = 1,2 



r\(V + rA 3 ) - n(V 3 + bA 3 ) i = 3 
Comparing Eqs. (|3.10|) and (|3.12j ), we see that 

J 3 + T 3 = 



(3.12) 



(3.13) 



which is a reflection of the axial symmetry of the solutions. 

Additional "body-frame" components can be constructed by utilizing contractions with 
the e£ ■ Components involving two such contractions vanish; i.e.,0 



On the other hand, the charges 



5 M _(») c _ n 



L ar — Pa e B °AB 



(3.14) 



(3.15) 



and their complex conjugates t ar are in general nonzero; they are related to the last term in 
Eq. (P3D by 



dw T dw 1 
~dt ~dl ~^f^ 3 



ti r ti r tir^lr 

b + r b — r 
Finally, we define a conserved relative U(l) charge 



(3.16) 



1 



1 



Q dt = -j= S A a dt = ~^ 
= (b\-K)(V + rA 3 



^AB + ^AB 



Sab dt= —7=^2 Pb ) Pa ) £ ab dt 
V2 a=1 



(3.17) 



Using Eqs. ( |2.23| ), ( j3.10|) , ( |3.12| ), ( |3.16| ), and (|3.17| ), we can now rewrite Q ang in terms of 
the body-frame components, obtaining 



1 



E 



bJl + lb + ^- J if - 2Vb^^J t T l 

N 



\t 2 {k + fib) fr{ 
+ ^b)[« + 2 ) + 2r/ 3 Q] + ^ + ii: 



r=3 



t\ T t\ r t2 r t2r 



b + r b — r 



dt 2 . (3.18) 



8 This is just a restatement of the fact that (N — 2) 2 eigenvalues of the moment of inertia tensor 
vanish because of invariance under a U(iV — 2) subgroup. 
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It must be remembered that the "body frame" components that appear in this expression 
are not in general conserved. While the square of the angular momentum, 



(3.19) 



the SU(iV) Casimir, 



N 



T = T{ + T 2 + T 2 + ^ (tlrtir + t2 r t2r) 



(3.20) 



r=3 



and the U(l) charge Q are constant, determining the time-dependence of the individual 
components of J and T is in general rather complex. There is considerable simplification, 
however, if the SU(iV) charges all vanish. The only nonzero charges remaining are then J 1; 
J 2 , and Q, with the former two entering £/ ang only through J 2 = J 2 + J|; using the relation 
between Vi and A4 that follow from the vanishing of the Tj, one finds that this is given by 



K 



+ r\6 2 + sin 2 6, 



(3.21) 



Equation (|3.18|) then reduces to 



Gs.ng(T = 0) 



6J 2 



t 2 (k + pib) 



+ 



a b 
— H 

k 2 K(b 2 — r 2 ) 



Q 2 dt 2 



(3.22) 



IV. MONOPOLE TRAJECTORIES 

We are now ready to obtain the equations of motion of our system and study the behavior 
of the cloud and the massive monopoles in scattering processes. We restrict ourselves here 
to the case of vanishing SU(iV) charges, T = 0, leaving a discussion of the more general case 
to the Appendix. 

Because the massive monopoles have well-defined positions while the massless ones do 
not, there may be a tendency to view the various scattering processes as resulting from inter- 
actions between the massive monopoles, with the cloud dynamics being in the background. 
This would be incorrect. The massive monopoles correspond to two mutually commuting 
SU(2) subgroups of SU(iV + 2). Were it not for the presence of the cloud, they would not 
interact at all and would move on straight lines with constant velocity. Hence, when follow- 
ing the motion of the massive monopoles it is important to also keep track of the size and 
motion of the cloud at the same time. 

Equations ( |2.22| ) and fl3.22|) provide a suitable starting point for obtaining the equations 
of motion of our system. Because the angle and phase variables enter the latter only through 
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conserved quantities, we can view Gang as defining an effective potential for r and b that can 
be combined with the first two terms in Eq. ( |2.22j ) to yield an effective LagrangianP] 



^eff 



Vf 2 + - 

2 4 



b + r 



6J 5 



bQ 2 



2r 2 (K + jjb) 2n(b 2 -r 2 ) 



(4.1) 



We begin our analysis of the solutions of this Lagrangian by examining the asymptotic 
solutions at large r and b; because the terms containing J and Q both fall rapidly in this 
regime, it is sufficient to do this asymptotic analysis with J = Q = 0. We find that the 
cloud and the massive monopoles decouple from each other at large times. After studying 
the properties of these asymptotic solutions, we then turn to the behavior at finite time, 
where the actual interactions take place, beginning with J = Q = 0. Finally, we consider 
the effects of nonzero J and Q. 

Before beginning this analysis, we must address the apparent singularity in the kinetic 
energy at r = b. This singularity is not physical, and can be eliminated by defining 



x 



s/k yb + r + y/b — r 
\Jb — r 



< y < x . 



(4.2) 



y = \fH y/b + r 

In terms of these coordinates, the Lagrangian takes the nonsingular form 



{x 2 + y 2 ) + 



2o? 



(xy + yx) 2 - 



a 2 (x 2 + y 2 ) 



where 



j2 2(x 2 + y 2 ) 
2x 2 y 2 (a 2 + x 2 + y 2 ) (x 2 — y 2 ) 2 

2k 



Q 2 



(4.3) 



(4.4) 



As we have indicated in Eq. (|4.2j ), the entire physical range 0<r<6<oois mapped onto 
the octant < y < x of the x-y plane. Nevertheless, it will become clear that if either J or 
Q vanish some trajectories can cross the boundaries of this octant. These can be understood 
as follows. The x-axis corresponds to r = 0. A trajectory crossing this line corresponds to 
a motion in which the two massive monopoles approach head-on, meet (at y — 0), and then 
pass through one another. To accommodate this, we extend the definition of x and y into 
the next octant by 



x 



y/n \Jb — r + y/b + r 



Vb 



VbT 



< — y < x . 



(4.5) 



9 We omit the term (/j,/k 2 )Q 2 , which gives a Q-dependent increase in energy but has no effect on 
the evolution of b or r. 
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Similarly, a trajectory crossing the line x = y corresponds to a motion in which the cloud 
shrinks to its minimum size, b = r, and then begins to grow again. To describe this, we 
define 



X = yfk~ y/b + r — yfb 



y 



y/K Vb + r + VF 



< x < y . 



(4.6) 



Extending these definitions in a similar fashion to the other octants, we obtain an eight-fold 
mapping of the physical r-b space onto the x-y plane. Throughout the plane we have the 
relations 



1 

2k 
1 



xy\ 



b=-(x 2 + y 2 ) 
4k 



(4.7) 



so that the x- and y-axes correspond to r = while the lines y = ±x correspond to minimum 
cloud size, b = r. 



A. Asymptotic solutions for J = Q = 



If J = Q = the Lagrangian of Eq. (|4.3|) reduces to 



L = ^{x 2 + y 2 ) + ^{xy + yx) 2 . (4.8) 
The resulting Euler-Lagrange equations are 

= (a 2 + y 2 )x + xyy + 2yxy 

= (a 2 + x 2 )y + xy'x + 2xxy . (4.9) 

Multiplying the first of these by x and the second by y, and then subtracting, we find that 
xx = yy. Using this to rewrite Eq. (|4.9|), we obtain 



= (a 2 + x 2 + y 2 )x + 2yxy 

= (a 2 + x 2 + y 2 )y + 2xxy . (4.10) 

There are two constants of the motion. The time-independence of the Lagrangian implies 
the conservation of the energy (which is actually equal to the Lagrangian) 

E = ^(x 2 + y 2 ) + ^(xy + yx) 2 . (4.11) 
In addition, multiplication of the first of Eq. fl4.10|) by y and the second by x shows that 
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B = ^xy(a 2 + x 2 + y 2 ) (4.12) 
a 1 

is also conserved. (We have no deep understanding as to why B is constant; we have not 
found any generalization of B that is conserved when either J or Q is nonzero.) 

If B ^ 0, neither x nor y can vanish, implying that both x and y vary monotonically; if 
B — but E ^ 0, then one of x or ?/ vanishes for all t, while the other never vanishes. Hence, 
there are no closed orbits in the x-y plane. It is then straightforward to show that x 2 + y 2 
must tend to infinity as t — > ±oo. There are two possibilities for the large time behavior. 
The first is that dy/dx tends to a nonzero finite constant. Examination of Eq. ( |4.10| ) shows 
that this is possible only if dy/dx — > ±1, which can only happen if E = ±B. This gives 
solutions where x 2 = y 2 varies linearly with time, corresponding to a solution with a minimal 
{b = r) cloud surrounding two massive monopoles with constant relative velocity. 

The other possibility is that \dy/dx\ tends to either or oo; because of the eight-fold 
mapping of b and r onto the x-y plane, the two are are equivalent. We will analyze in detail 
the case where dy/dx — > 0, so that asymptotically \x\ ^> \y\. By integrating the dominant 
terms in the second of Eq. ( |4.10| ), we find that yx 2 asymptotically tends to a constant. 
Substituting this result into Eq. ( |4.12| ) for B, we see that x must tend to a constant. It is 



then straightforward to obtain the asymptotic solution 

2,;v °it-t ) + o(r') 



x 



yo 



y = y -^ + O(r 2 ) (4.13) 

where yo, t , f , and (3 are constants of integration. Of these constants, t simply corresponds 
to a shift of the zero of time, while (3 is related to the others by 

y^ + yl)B 

H 4kv E v ' 

The meanings of Vq and yo are clarified by writing the asymptotic solution in terms of b 
and r: 

r = v \t\ + . . . 

b =^t 2 + .... (4.15) 
Vo 

We see that asymptotically the massive monopoles move with a constant relative velocity 
v , while b varies quadratically with time.Q Substituting this asymptotic solution back into 
the energy, we find that 



10 This would imply that as t approaches oo (— oo), the cloud parameter increases (decreases) at 
a rate that exceeds the speed of light; we will return to this point below. 
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+ ^ + (4-16) 

As t — > ±00 the terms indicated by dots vanish, while the two terms shown explicitly tend 
to constants. These latter two terms can be interpreted as a massive monopole energy E r 
and a cloud energy E b , with 

V 2 o = a 2 ^. (4.17) 

Because of the multiple mapping of b and r into the x-y plane, for each asymptotic 
solution in Eq. (|4.13| ) there are seven other physically equivalent solutions that are obtained 
by reversing the signs of x or y or by the interchange of x and y. 

B. Joining the asymptotic solutions when J = Q = 

Having found the possible asymptotic behaviors, we now want to obtain solutions of 
Eq. ( |4.10| ) for finite times. As with the asymptotic solutions, these will depend on four con- 
stants of integration. One of these can be absorbed in a rescaling of t (with a corresponding 
rescaling of E and B) and a second in a shift of the zero of time. Hence, the trajectories 
in the x-y plane depend nontrivially only on two constants, which we take to be B/E and 
E r /Eb] without loss of generality, we can choose the trajectories to begin in the first octant, 
so that the latter quantity is equal to y 2 /a 2 . Equations ( 4.11|) and ([4.12 ) can be combined 



to yield a quadratic equation for dy/dx = y/x. The solution of this equation gives dy/dx 
as a function of x, y, and B/E. Hence, two trajectories with the same value of B/E but 
different yo cannot cross. While our numerical solutions indicate that the same may be true 
for trajectories with the same yo and different B/E, we have no proof of this statement. 

We used numerical integration to obtain solutions for various values of yo and B/E. We 
will begin by describing our results for y = a (i.e., initially equal values of E r and Eb). In 
Fig. 1, we show some trajectories in the x-y plane for y = a and several values of B/E. 
Figure 2 shows r and b as functions of time for these trajectories. Note that r and b have 
been plotted in units of r c = 2k/ fi; recall that this is approximately equal to the sum of the 
core radii of the two massive monopoles. 

In examining these solutions, it is useful to focus on the points where b = r (corresponding 
to crossing one of the 45-degree lines x = ±y) and where r = (corresponding to crossing 
either of the coordinate axes in the x-y plane). When b = r, the cloud is of minimal sizej^j 



11 This is not necessarily the minimum value of b. In most cases r < b when b reaches its minimum, 
so that the smallest value of b corresponds to a nonminimal cloud. 
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-5 5 -1 1 

x/a x/a 
FIG. 1. Several trajectories in the x-y plane for yo = a. Beginning from the top, these 

correspond to values of B/E of -1.40, -0.60, 0.66, 0.85, 1.20, and 3.00. The graph on the right is a 
blow-up of the region near the origin. 

and can be viewed as having shrunk to the line segment joining the massive monopoles. The 
field configuration is an embedding of a purely Abelian SU(3) — > U(l) x U(l) solution, and 
the component of the field in the unbroken SU(iV) is a pure dipole field. In terms of massless 
monopole positions, this corresponds to a configuration where all of the massless monopoles 
are located on top of one of the massive monopoles. All trajectories pass through at least 
one such minimal cloud configuration. If \B / E\ < 1, the trajectory passes through two such 
configurations. Whether b = r once or twice does not appear to have a dramatic influence 
on the large time behavior of the solutions. 

The points on a trajectory where r = are, of course, the points where the two massive 
monopoles pass through one another. Since x and y both vary monotonically, it is clear 
that every trajectory must cross at least one coordinate axis and so must have at least one 
such point. Thus, the two massive monopoles that are approaching each other on a straight 
line (because J = 0) cannot stop and reverse direction before they meet. However, it is 
possible for them, after passing through each other once, to stop, reverse direction, pass 
through each other a second time, and finally emerge with their original direction of motion 
reversed. Whether this second possibility is realized is determined, for a given y , by B/E. 
For y = a, the monopoles reverse direction if 0.51 < B/E < 2.60; the boundaries of this 
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FIG. 2. Plots of b (dashed line) and r (solid line) as functions of t for the trajectories shown 
in Fig. 1. The progression from (a) to (f) corresponds to moving from top to bottom in Fig. 1. 
Time is shown in units of a/2y/~E, while b and r are plotted in units of r c = 2k/ fi. Note that the 
scale varies from plot to plot. 
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FIG. 3. Energy transfer from the cloud to the massive monopoles as a function of B/E, for 
Vo = a. 

range correspond to solutions where the massive monopole velocities vanish as t — > oo. 

The trajectories displayed in Figs. 1 and 2 were chosen to illustrate the various possible 
behaviors (b — r once or twice, r = once or twice). Several points should be noted. First, 
there does not appear to be any particular striking effect when r = 0. Because r is by 
definition positive, f changes sign, but the magnitude of f varies continuously. Similarly, 
b(t) shows no unusual behavior when r vanishes. All this is consistent with our expectations; 
since the two massive monopoles have no direct interactions, they should be able to pass 
through each other with little effect. 

By contrast, each of the massive monopoles does interact directly with one of the massless 
monopoles comprising the cloud. This might suggest that there should be a noticeable effect 
when a massive and a massless monopole position coincide; i.e., when b = r. On the other 
hand, the massless monopoles are hardly point particles, since the entire cloud region can 
be viewed as being composed of the massless cores. Nevertheless, it does seem that the 
greatest effect on the trajectories occurs when b r. Thus, the largest changes in the slope 
of r(t) coincide almost exactly with the vanishing of b — r. Note also that the effect of these 
massive-massless "collisions" seems to be greatest when they occur at small b and r. This 
is perhaps clearest in the trajectories where the interaction is strong enough to reverse the 
direction of motion of the massive monopoles, so that r = twice. In these cases, b and r 
coincide when r is so small that the two massive cores overlap each other. 

Even when the massive monopoles do not change direction, we can gauge the strength of 
their interaction with the cloud by the amount of energy transferred between the cloud and 
the massive monopoles. In Fig. 3 we show this as a function of B/E for y = a. Although 
the details of this plot vary somewhat with y , two points are fixed: 

1) If B/E — 0, y vanishes identically, so y is constant in time and there is no energy 
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-4-2 2 4 

B/E 

FIG. 4. Minimum cloud parameter 6 m i n (dashed line) and the corresponding value of 6 m j n — r 
(solid line) as functions of B/E, for yo = a - Both b and r are plotted in units of 2k/ fi. 

transfer. 

2) If _B/_E = ±1, the trajectory asymptotically approaches one of the lines x = ±y at 
large time, corresponding to a final state with minimal cloud, b = r. In this case, all of the 
initial cloud energy is transferred to the massive monopole kinetic energy. 

Figure 3 should be compared with Fig. 4, where we we plot the minimum value of b as a 
function of B/E. We see that the largest energy transfers are associated with small values 
of b, and that the energy transfer decreases as b min increases for \B/E\ ^> I. In Fig. 4 we 
have also shown the value of b — r at the time that b achieves its minimum. Note that this 
value is always less than or comparable to the monopole core radius. 

We have also studied solutions with y = 100a (i.e., E r 3> E b ) and y = 0.01a (E r <C E b ). 
The general picture is very much as for the case of y = a. The analogues of Figs. 1 and 2 
are qualitatively quite similar. We again find that b — r is small when the cloud parameter 
b achieves its minimum size, and the greatest amount of energy transfer between the cloud 
and the massive monopoles occurs when 6 min is small. One notable difference is found in 
those yo = 0.01a solutions where the massive monopoles reverse direction. In contrast with 
the yo = a case, the massive monopoles need not be overlapping at the time of this reversal. 
This is readily understood. Because the cloud energy is so much greater than the massive 
monopole energy, relatively little energy transfer is required to reverse the sign of f, and so 
b at the time of the collision can be larger than it was in the previous case. This emphasizes 
that the critical factor for the yo = a reversal was that b was small; the fact that the massive 
cores overlapped was coincidental. 
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x/a x/a 
FIG. 5. Several trajectories in the x-y plane for nonzero J. The values of yo and B/E are the 

same as in Fig. 1, while J = 0.001a^/E. The trajectories in the left box correspond to the top three 
trajectories in Fig. 1, while those in the right box correspond to the bottom three trajectories in 
Fig. 1. 

C. Solutions with nonzero J and Q 



The effective potential terms in Eq. ( j4.1|) come into play when either J or Q is nonzero. 
With nonvanishing angular momentum, there is a centrifugal barrier that prevents the van- 
ishing of r. The trajectory in the x-y plane cannot cross either axis, and so is confined to a 
single quadrant. For large J, this barrier ensures that r, and hence b, is always large, and 
thus suppresses interactions between the massive monopoles and the cloud. For small but 
nonzero J, the effect of the centrifugal barrier is to make the trajectories appear to reflect 
off the x- and y-axes. We illustrate this in Fig. 5. Except for the value of J, the initial data 



for the trajectories in this figure are the same as for those shown in Fig. 1. Equation (|3.19| ) 
for J 2 implies that the motion of the massive monopoles is confined to a plane. In Fig. 6, 
we illustrate this motion by showing the path of r in the plane perpendicular to J for one of 
the trajectories shown in Fig. 5. Note once again that there is little evidence of interaction 
when the massive monopoles are at their closest approach. 

If Q ^ 0, there is a barrier that forbids the vanishing of b — r, or equivalently, the 
crossing of the lines x = ±y. When Q is large, b — r can never be small, and so there is little 
interaction among the monopoles. For small nonzero Q, we find behavior analogous to that 
for J ^ 0, with trajectories appearing to bounce off the lines x = ±y. 

V. CONCLUDING REMARKS 

In this paper we have studied the properties of the magnetically charged counterparts of 
the electrically-charged massless particles that arise when a gauge theory is spontaneously 
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FIG. 6. A portion of the orbit of the massive monopole separation r in the plane perpendicular 
to J. Note that the turnaround of the monopole motion does not occur at the time of closest 
approach. The trajectory shown here corresponds to B/E = 1.2, yo = a, and J = 0.002aV^E. 



broken to a non-Abelian subgroup. Previous studies of static BPS solutions have shown 
that these massless monopoles coalesce into a cloud of non-Abelian field surrounding one 
or more massive monopoles. In the scattering solutions that we have studied, these clouds 
act very much like thin shells, with a size of order b, that contract until they collide with 
and bounce off the cores of the massive monopoles. At large times, when this shell is far 
from the massive cores, the cloud and the massive monopoles are essentially decoupled, and 
have separately conserved energies. Their interaction occurs primarily at the time that the 
shell overlaps the massive cores. We find that the strength of this interaction, as measured, 
e.g., by the amount of energy transfer, tends to be greatest if b is small at the time of this 
overlap. 

Further evidence in support of this picture can be obtained from another system that 
has been studied in some detail. When SU(3) gauge theory is spontaneously broken to 
SU(2)xU(l) there are two species of fundamental monopoles, one massive and one massless. 
The Nahm data and moduli space metric for the BPS solutions containing one massless and 



two massive monopoles were obtained by Dancer [|14| , |15|| . These solutions depend on twelve 
parameters. Ten of these correspond to the center-of-mass position, global SU(2)xU(l) 
phases, and overall spatial rotations. The remaining two, k and D, arise as parameters in 
elliptic functions that enter the Nahm data. These span a geodesic submanifold which, after 
a nonlinear change of variables to eliminate coordinate singularities in the moduli space 
metric, is often illustrated as in Fig. 7. (This figure is actually a six-fold covering of the 
submanifold, completely analogous to our eight-fold mapping of b and r onto the x-y plane.) 

The precise relations between position in this figure and cloud size and monopole separa- 
tion are rather complex, with analytic expressions only known asymptotically . However, 
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FIG. 7. The two-dimensional Dancer submanifold discussed in the text. The hyperbolic solu- 
tions correspond to points on the solid lines, while the trigonometric solutions lie on the dashed 
lines. The boundaries, indicated by dotted lines, correspond to embedded SU(2) two-monopole 
solutions. 



the qualitative picture is easily described. Points far out along the three legs of the diagram 
correspond to solutions containing two well-separated massive monopoles, with the separa- 
tion increasing with distance along the leg, while points in the central part of the figure 
correspond to solutions in which the massive monopole cores overlap and deform each other. 
The boundaries of the figure correspond to embedded SU(2) two-monopole solutions; in our 
language, these are solutions with infinite cloud parameter. (These boundaries are infinitely 
far, in metric distance, from any interior point.) The lines that bisect the legs correspond 



to special axially symmetric solutions fll5 |. The solid portions of these lines in Fig. 7 cor- 
respond to "hyperbolic solutions" that can be viewed as solutions with minimal cloud (the 
analogue of our b = r solutions) and with the separation between the massive monopoles 
increasing with distance from the center. The dashed portions of these lines correspond to 
"trigonometric solutions" that are composed of two overlapping massive monopoles and a 
massless cloud that increases from minimal size at the center to infinite size at the boundary 
of the figure. 

Geodesies on this two-dimensional submanifold correspond to scattering solutions with 
vanishing angular momentum and internal charges. In contrast to the SU(iV) example we 
have studied in this paper, there are direct interactions between the two massive monopoles. 
Consequently, there is nontrivial scattering even when the cloud is infinitely large. Thus, 
each of the three geodesies that bound the moduli space in Fig. 7 corresponds to a process 
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in which the two massive monopoles start infinitely far apart, approach head-on, and then 
recede at right angles to their initial motion; during the entire process the cloud remains 
infinitely large. 

Numerical solutions of the geodesic equations were obtained by Dancer and Leese fllTf . 
Examining these, we see behavior quite consistent with our interpretation of our SU(iV) 
solutions, in that the interaction between the cloud and the massive monopoles is greater if 
the massive monopoles are closer together at the time of the minimal cloud configuration. 
We see this in two ways. First, there are some geodesies that start in one leg of the figure, 
cross the line of trigonometric solutions and begin to move down another leg, and then, 
after crossing the line of hyperbolic solutions, cross a second line of trigonometric solutions 
and exit via third leg. These may be viewed as analogous to the SU(iV) scattering solutions 
in which the two massive monopoles reverse direction. As with the latter solutions, these 
"three-leg trajectories" are found only if the massive monopoles are sufficiently close in the 
minimal cloud configuration; in other words, if the trajectory goes too far down the second 
leg, the massive monopoles cannot reverse direction. 

The effects of cloud interactions on the "two-leg trajectories" that cross a line of hyper- 
bolic solutions but only a single line of trigonometric solutions are more subtle. At large 
times these trajectories lie between the boundary and the line of hyperbolic solutions, with 
the cloud size increasing as the relative distance to the boundary decreases. Comparing 
the distance of the geodesies from the boundary at points equally far down the initial and 
final legs thus gives a measure of the net interaction between the cloud and the massless 
monopoles, somewhat analogous to the SU(iV) energy transfer that we plotted in Fig. 3. 
Examination of the Dancer-Leese geodesies shows that the greatest change between initial 
and final legs, and thus the greatest interaction, occurs if the geodesic crosses the line of hy- 
perbolic solutions nearest the center of the figure; i.e., nearest the point of minimal massless 
monopole separation. 

Before closing, we must address the validity of the moduli space approximation. There 
is a potential problem in this regard in any theory where the excitation of massless particle 
modes is a possibility. The fact that the asymptotic solution of Eq. (|4.15|) has the cloud 
parameter b increasing faster than the speed of light at large times is an indication that there 
actually is a breakdown of the moduli space approximation in the example we have studied. 
Considerable insight in this issue can be obtained by studying [Tj| an SO (5) example with 
a spherically symmetric cloud. This symmetry makes it possible to compare the predictions 
of the moduli space approximation with a numerical solution of the full field equations, 
starting from a slowly varying initial configuration. It turns out that the moduli space 
approximation continues to be reliable until b becomes of order unity. After that time, the 
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edge of the cloud develops into a rather well-defined front expanding outward at the speed 
of light. Behind this front, in the cloud interior, the time evolution of the fields continues 
to be well-described by the moduli space solution (i.e., as if b ~ t 2 ). Ahead of the front, the 
evolution of the fields is much slower, consistent with a slow linear growth of the cloud at a 
rate fixed by the initial value of b. 

How does this affect the picture of monopole-cloud scattering that we have developed? 
Two cases can be distinguished. In one, a static solution with finite r and b is subject to 
an external disturbance that induces small, but nonzero, values for r and b. As long as b 
remains small (energy conservation ensures that r will remain small) the moduli space de- 
scription of the interaction between the cloud and the massive monopoles should be reliable. 
Eventually b will become of order unity, but since r ^ b when this happens, the moduli 
space approximation should continue to give a good description of the subsequent motion 
of the massive particles. Thus the main modification to our picture is that the velocity of 
the cloud asymptotically approaches unity, rather than increasing without bound. 

In the other case, the initial low-energy configuration approximates a solution with large 
cloud that is contracting with the magnitude of b being of order unity. We expect that the 
deviation from the moduli space solution will lead to the emission of long wavelength, low 
amplitude waves of the massless boson fields that propagate at the speed of light. This will 
leave behind a configuration whose evolution then proceeds much as in the previous case. 

In closing, we note that while our focus in this paper has been on the massless monopole 
clouds in classical multimonopole configurations, much of the underlying motivation was to 
gain insight into the states of the electrically-charged elementary particle sector that are 
dual to these. Thus, it would be of great interest to study the BPS and near-BPS states 
containing massless and massive elementary particles with non-Abelian charge with the aim 
of identifying the counterparts of the cloud parameter b or of the asymptotic separation of 
the energy between the cloud and the massless monopoles. 
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APPENDIX A: 

In this appendix we show explicitly how to express the metric Q rc \ of Eq. Q2.7| ) in terms 
of the Euler angles. We restrict ourselves here to the case of SU(4) — > U(l) x SU(2) x U(l). 
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A natural set of parameters is the cloud size b; the massive monopole separation r; three 
spatial Euler angles 9, <fi, and ip] three internal SU(2) Euler angles a, (3, and 7; and the 
relative U(l) phase %• This is a total of nine but, as we discussed in Sec. 2, we only expect 
eight parameters. We will see later that this conflict is resolved by the identification of a 
spatial Euler angle with one of the internal SU(2) Euler angles. 

The z a A transform as doublets under both the spatial and the internal SU(2), with the 
upper and lower indices being spatial and internal, respectively. Hence the z A for arbitrary 
values of the spatial Euler angles must be of the form 



z 




U{0,M)[, A \ (Al) 



where a A and b A are functions of b, r, Xi an d the internal Euler angles. We choose a 
parameterization of U (8, <p, ip) such that 

z a A = a A v a + b A v a (A2) 

where 

'cosfeiH^) 
sinfesHH-*) 

/ sir, iMih-*) \ 
«»= * ft i!t ^ (A3) 

Similarly, a A and b A can be written as linear combinations of 



u A 



sinfet(-^ +/3 ) 
cos ^e^~ 7_/3 ^ 



with coefficients that depend only on r, b, and \. By performing a gauge rotation, we can 
make a A proportional to u A . We then have 

z a A = Au A v a + Bu A v A + Cu A v A (A5) 

The 2 relative phases among A, B and C can be absorbed in redefinitions of 7 and ip. 
Then, by a redefinition A — > Ae^ x , B — > Be% x and C — > Ce^ x , we can make the A, B, C all 
real and non-negative. Equations ( |2.5| ) and ( |2.6| ) then give the constraints 

A 2 + C 2 = ^(b + r) 

B 2 = \{b-r) 

BC = (A6) 
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These are solved by A = \\Jb + r, B = -yb — r, C = 0. This gives us 



z, = W sm — cos -e w + W cos — sin-e 2iWJ e2^ +x w , 



, " L/WQ — Dill ~G" 

2 2 2 J 



= n/ cos — cos -e zvy+w _ , / sm _ sm _ e2 a+w e2 >. 



sm — sm-e _ . cos _ cos _ e2 m-w ga^t-x-t-w 



2 2 2 



b _±L cos I sin ?e-J(TH) + , /^I s in - cos -e^+<« I c *H»+**>. (A7) 
2 2 2 y 2 2 2 / 

Notice that 7 and ^ always come together in the combination 7 + ■?/>. This is a manifestation 
of the axial symmetry of this configuration: for any spatial axial rotation that changes the 
parameter x/j, one can do an opposite axial gauge rotation that changes 7, so that the metric 
is invariant. 

Substituting these relations into the Eq. Q2.71) gives 



'rel 



where 



+ -(/ir 2 + Kb)(al + a 2 2 ) + -/t6(r 2 + r|) + «V& 2 - r^-o^n + a 2 r 2 ) 
+ 2{ b -^)^ + T ^-^ ( " 3 + T ^ + 2(^5)* 2 < A8) 



o"i = — sin ipdO + cos ^ sin 6 1 cf</>, 
o"2 = cos ip d9 + sin -0 sin 6* d0, 

cr 3 = d0 + cos # (A9) 



and 



ti = — sin ■yda + cos 7 sin a d/3, 
T2 = cos 7 da + sin 7 sin a d{3, 

r 3 = c/7 + cos a d/3, (A10) 

are the left-invariant one-forms for the rotational and internal SU(2)'s, respectively. 

There are several interesting limits of this metric. If the relative distance between the 
two massive monopoles r and the relative U(l) charge Q both vanish, it reduces to the 
metric for one massive and one massless monopole in SO(5) |^,|12|: 



28 



'SO(5) 



K 

2 



+ b{rl + r 2 2 + rl) 



(All) 



In the limit where we have minimal cloud b = r and vanishing internal SU(2) charges, 



we get the metric for two distinct massive SU(3) monopoles [nj|. The angles ip and 7 can 
be absorbed into a redefinition of %■> an d we obtain 



1 



£SU(3) = o Z 1 + - rfr + r (°1 + °" 



2 k + jir 



-dx + cos 6 dcj)) . 



(A12) 



We can also consider the limit where fi is much greater than both n/b and the initial 
energy of the cloud. To leading order, the massive monopole motion decouples from that of 
the cloud, and we may set r = 0. Making the change of variables b = \ and r = f, we get 



the Eguchi-Hanson |T{| metric 12 . 



/EH 




2 1 2 



1 



4C 



1 'J 



(A13) 



When C = 0, this reduces to Eq. ( |A11| ). 

In order to study the classical dynamics for the general case with nonzero angular mo- 
menta and internal charges, it is most convenient to first eliminate the redundant parameters. 
This can be done, e.g., by simply choosing %p = 0. We then write the angular part of the 
metric as 



?ang 



Iij(r, fyufu 3 



(A14) 



where the Ui (i = 1, . . . , 6) are one-forms expressed in terms of the various rotational and 
gauge angles. These satisfy an algebra of the form 



2 3 



(A15) 



where the Cj k are antisymmetric in j and k. We then define the charges to be Xj = I^uj- K 
Using Eq. (|A15 ), we can write the angular equations of motion for Eq. ( |A14| ) as 

d 



di x t + ci l (i- 1 ) k %x l = o. 



(A16) 



Ordinarily, we could have chosen the u J ' to include the three a a and the three r a . However, 
after eliminating the redundant angular variable, this is not possible. Instead, we choose 



CO 



1 - - sin 6 , 



00 



2 - d9. 



LO = Ti, 



12 



We thank Kimyeong Lee for pointing this out to us. 
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u A = r 2 , u b = r 3 + cos 9 d<j>, uj b = dx ■ (A17) 

With this choice, the Cj fc are not constants, as they are in more familiar examples. The 
nonzero Cj fc are given by 

C21 ~~ ^41 = ^13 = Cot 

C| 5 = C53 = C\ x = C| 4 = 1 . (A18) 



Iij can be obtained from Eq. ( |A8|) for the metric. It is a block diagonal matrix formed by 
three two-by-two blocks, with the first two blocks being identical. It is straightforward to 
verify that Eq. (|A16|) implies that the square of the angular momentum, Eq. ( |3.19| ), and the 
square of the internal SU(2) charge, Eq. (|3.20| ), are indeed constant. 
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